Factorization of shell-model ground-states 
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We present a new method that accurately approximates the shell-model ground-state by products 
of suitable states. The optimal factors are determined by a variational principle and result from the 
solution of rather low-dimensional eigenvalue problems. The power of this method is demonstrated 
by computations of ground-states and low-lying excitations in sd-shell and p/-shell nuclei. 
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Realistic nuclear structure models are notoriously diffi- 
cult to solve due to the complexity of the nucleon-nucleon 
interaction and the sheer size of the model spaces. Exact 
diagonalizations are now possible forp /-shell nuclei 0,0 
and for sufficiently light systems 0, Q, and Quantum 
Monte Carlo calculations [a, |(| have solved light nuclei 
up to mass number A — 10. For cases where an exact di- 
agonalization is not feasible, various accurate approxima- 
tions are employed. Important examples are stochastic 
methods like Shell-Model Monte Carlo @, H and Monte 
Carlo Shell-Model 0. Recently, non-stochastic approx- 
imations have been developed. Examples are the Ex- 
ponential Convergence Method d |H El E3 and the 
application of density matrix renormalization group to 
nuclear structure problems 0, 0, 0. These latter 
methods truncate the Hilbert space to those states that 
accurately approximate the ground-state, and the results 
of such calculations usually converge exponentially as one 
increases the number of retained states. The correct 
identification of the most suitable states clearly becomes 
crucial for these approximations to be useful. In this 
work, we propose a solution to this problem and present 
a method that obtains the optimal states from a varia- 
tional principle. 

We divide the set of single-particle orbitals into two 
subsets, V and A/", respectively. Here V represents the 
set of all proton orbitals while M denotes the set of all 
neutron orbitals. We label many-body basis states within 
each subset as {|7r a ),a: = l,...,dp} and {|^ Q ),cv = 
1, . . . , d_\f}. We expand a nuclear many-body state 
in terms of the proton and neutron states 



(i) 



a 



This expansion is not unique since the amplitudes $ a p 
depend on the choice of basis states within the two sub- 
sets. There is, however, a preferred basis in which the 
amplitudes ^ a f3 are diagonal. This basis is obtained from 
a singular value decomposition of the amplitude matrix 
Vl/a/j, i.e. \& = USV T , with diagonal S and orthogonal 
matrices U and V. 
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FIG. 1: Singular values s 2 for ground-states of 20 Ne (*, d-p = 
66), 22 Ne (x, d v = 66), 44 Ti (+, d v = 190), and 24 Mg (•, 
d v = 495). 



Here Sj denote the singular values while the correlated 
P-states \pj) and the jV-states \hj) are orthonormal sets 
of left and right singular vectors, respectively. The non- 
negative singular values s\ > S2 > S3 . . . fulfill J2j s j = 1 
due to wave function normalization. It is interesting to 
take ground-states of realistic nuclear many-body Hamil- 
tonians and compute their singular value decomposition 
in terms of proton and neutron states. To this purpose 
we take the amplitude matrix ^ a p of the ground state 
and perform a numerical singular value decomposition of 
this matrix. Figure ^ shows the squares of the singular 
values for the ground-states of sd-shell nuclei 20 Ne, 22 Ne 
and 24 Mg (from the USD interaction ^tJ) and for the pf- 
shell nucleus 44 Ti (from the KB3 interaction [lj|)- Note 
that the singular values decrease rapidly (degeneracies 
are due to spin/isospin symmetry). This suggests that a 
truncation of the factorization (J2J should yield accurate 
approximations to the ground-state. The density matrix 
renormalization group (DMRG) [l9|, exploits this rapid 
fall-off of singular values in a wave-function factorization. 
For first applications of this method to nuclear structure 
problems we refer the reader to Refs. 0, 0, 0] . 
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In this paper, we present an alternative technique that 
efficiently factorizes shell-model ground-states and deter- 
mines the optimal factors for a given truncation. We 
make the ansatz 



(3) 



3=1 



Here, the unknown factors are the P-states \pj) and the 
A/"-states \rij). The truncation is controlled by the pa- 
rameter 51 which counts the number of desired factors. 
Figure ^ suggests that 51 <C min (dp, d_\f) yields accurate 
approximations to shell- model ground-states. This is also 
the result of our numerical computations below. 

Let H be the nuclear many-body Hamiltonian. The 
unknown P-states \pj) and TV-states \rij) in Eq. © 
are obtained from a variation of the energy E — 
(V>|#|V>/(^#), which yields (j = 1, . . . , 51) 

Q 

^({nAH^-E^mj) \ Pi ) = 0, 

i=l 
Q 

Y,(<Pi\H\Pi)-E<pj\pi)) \m) = 0. (4) 

i=l 

These nonlinear equations are not easy to solve simul- 
taneously. Note however that for fixed neutron (proton) 
states the first (second) set of these equations constitutes 
a generalized eigenvalue problem for the proton (neutron) 
states. Consider the first set of the Eqs. (0J. The opera- 
tor (m\H\nj) acts on P-space and is determined by the 
nuclear structure Hamiltonian 



H = H A r + H v + VpM, 
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with 
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Hp = ^£ p aja p + - X! v pp'q'q^W P '^q'^ ( 6 ) 
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VvH = ^2 v pnp'n'd\a { n a nl h vl . 
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Here, we use indices p, q and to, n to refer to proton and 
neutron orbitals, respectively. The antisymmetric two- 
body matrix elements are denoted as Vijki- 

Thus, the P-space Hamilton operator (fij\H \rii) is 



(ni\H\nj) = ^ ^Vpnp'n'intldlan'lnj) d\a p > 

p,p' \n,n' / 

+ (ni\Hx\nj) + (m\rij)H-p. (7) 

Note that the neutron- proton interaction V-pjj results 
into a one-body proton operator while the neutron 



Hamiltonian Hj^ yields a constant. This concludes 
the detailed explanation of the first set of equations in 
Eq. J2J. The second set has an identical structure, only 
the role of neutrons and protons is reversed. 

We solve the coupled set of nonlinear equations Q 
iteratively as follows. We choose a random set of ini- 
tial TV-states {\rij),j — 1, . . . , 51} and solve the first set 
of Eq. igj for those P-states {\ P j),j — 1, . . . , 0} that 
yield the lowest energy E. We then use this solution for 
the second set of Eq. which yields improved TV-states 
{\rij},j = 1, . . • , 51}. We iterate this procedure until the 
energy E is converged, typically 5-20 times. The ad- 
vantage of the factorization method is that the dimen- 
sionality of the eigenvalue problem is 51 x max (d-p, dj\f) 
while an exact diagonalization scales like dp x dj\f. Be- 
fore we present numerical results we note two further 
developments of the method. First, we will exploit the 
non-uniqueness of the ansatz J3J to reduce the general- 
ized eigenvalue problem to a standard eigenvalue prob- 
lem. Second, we will use the rotational symmetry of the 
shell-model Hamiltonian (J^J to give a m-scheme formu- 
lation of the factorization method. 

To reduce the Eqs. JIJ to a standard eigenvalue prob- 
lem we choose an orthonormal set of initial TV-states 
{\rij),j = 1, . . . , 51} as input to the first eigenvalue prob- 
lem. This reduces the "overlap" matrix (rij\rii) to a unit 
matrix, and we solve a standard eigenvalue problem to 
obtain the P-states {\pj),j — 1,...,51}. The resulting 
P-states will not be orthogonal in general. Their coeffi- 
cient matrix C with elements C a j = {n a \pj} may, how- 
ever, be factorized in a singular value decomposition as 
C = UDV T . Here D denotes a diagonal 51 x 51 matrix 
while U is a (column) orthogonal dp x 51 matrix, and V 
is a orthogonal 51 x 51 matrix. The transformed states 



a=l 

K> = X>>*> 



i=l 

are orthonormal in P-space and TV-space, respectively, 
and fulfill 

\i') = J2\p^ = J2 D n\p'M)- (8) 

3=1 3 = 1 

We then input the new P-states {\p'j),j = 1,...,51} 
to the second set in Eq. 10}, which poses a standard 
eigenvalue problem for the TV-states. The transformed 
neutron states Djj\n'j) are good starting vectors for the 
Lanczos iteration. We iterate the whole procedure un- 
til the ground-state energy converges to an appropriate 
level of accuracy. Note that the singular value decomposi- 
tion is very inexpensive compared to the diagonalization. 
This procedure also has the advantage that it yields the 
singular values Djj of the ground-state factorization. 
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E [MeV] 
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E [MeV] 


d 


Factorization 


-10.74 


190 


-29.55 


4,845 


-50.41 


38,760 


-76.23 


125,970 


m-scheme Fact. 


-12.74 


190 


-31.06 


4,845 


-52.03 


38,760 


-77.1 


125,970 


Exact 


-13.88 


4,000 


-32.95 


1,963,461 


-54.27 


109,954,620 


-78.46 


1, 087,455,228 



TABLE I: 

Results for ground-state energies Eq of p/-shell nuclei obtained from the factorization with just one state and the most severe 
m-scheme factorization compared to the exact results from Ref. Q. The m-scheme dimension of the corresponding eigenvalue 

problem, d, is also listed. 



To include the axial symmetry ("m-scheme") into the 
factorization we consider only products of proton and 
neutron states that have angular momentum projection 
J z = 0. To this purpose we modify the ansatz © as 

w=£f>5 w0 >i»l _w0 >. ( 9 ) 

m j=l 

where |pj m) ), j = !,-■■, O m (|nj~ m) ), j = 1, . . . , Q m ) de- 
note many-proton states (many-neutron states) with an- 
gular momentum projection J z = m (J z = —to), and the 
sum over to runs over all value of J z that are realized. 
The ansatz ® again leads to a generalized eigenvalue 
problem similar to Eq. which we also cast into a 
standard eigenvalue problem by maintaining orthogonal- 
ity between P-states and between the TV-states through 
singular value decompositions. The number of factors 
used in the m-scheme factorization is given by the pa- 
rameters Q rn . In the following, we use 

fim(a) =max(l,aIC"), (10) 

where fi™ ax is the maximal dimension of the correspond- 
ing subspace. For a = the most severe truncation is 
obtained while a — 1 leads to an eigenvalue problem 
with the same dimension as an exact diagonalization in 
m-scheme. In the latter case, the first iteration (i.e. solv- 
ing the first set of Eq. (0} for the proton states and using 
random, orthonormal neutron states as input) leads to 
the exact solution. This is due to the non-uniqueness of 
the ansatz (0. Note that other abelian symmetries and 
point symmetries can implemented in a similar manner 
as the axial symmetry. This flexibility is important since 
it allows us to consider factorizations that differ from the 
neutron-proton factorization proposed here. In neutron 
rich nuclei, for instance, there might be a large imbal- 
ance between the sizes of the proton and neutron space. 
In such a situation one might consider a "mixed" factor- 
ization where part of the neutron orbitals and all proton 
orbitals constitute one factor space while the other fac- 
tor space consists of the remaining neutron orbitals. Note 
also that factorization method proposed here is not re- 
stricted to nuclear structure problems but could well be 
applied to problems in quantum chemistry or condensed 
matter. 



We turn to numerical tests of the factorization method. 
We solve the eigenvalue equations with the sparse ma- 
trix solver ARPACK [2jj. In a first step we set fi = 1 in 
the ansatz @ and compare the most severely truncated 
factorization with exact diagonalization. The factoriza- 
tion requires us to iteratively solve eigenvalue problems 
with rather low dimension d-p and dj\f, respectively. Ta- 
ble [I] shows the results for the p/-shell nuclei 44 Ti, 48 Cr, 
52 Fe, and 56 Ni. The results from the severely truncated 
factorization deviate about 2-4 MeV from the exact re- 
sults (from Ref. |2]). For comparison, we also list the 
dimension of the eigenvalue problem. It is also interest- 
ing to use the m-scheme factorization and set Q m — 1 in 
Eq. © . This corresponds to a = in Eq. I|1U[1 and yields 
a problem of equal dimension as the factorization (|3J) with 
truncation f2 = 1. Results are also listed in Table [I] The 
advantage of the m-scheme factorization is clearly visible 
as the corresponding results typically differ only 1-2 MeV 
from the exact results. This makes this severe truncation 
particularly interesting to compute good starting vectors 
for large-scale nuclear structure problems. After these 
encouraging results we will in the following include more 
factors to obtain more accurate ground-state factoriza- 
tions. 

We use the m-scheme factorization to compute accu- 
rate approximations to ground-state energies. As a first 
test case we consider the sd-shell nucleus 24 Mg with the 
USD interaction ^tJ- The exact diagonalization yields 
the ground-states energy E = —87.08 MeV and has m- 
scheme dimension d max = 28503. We solve the eigen- 
value problem |0J for the ground state but also record 
a few excited energy solutions. The hollow data points 
in Fig. |21 show the resulting excitation spectrum versus 
the dimension of the eigenvalue problem relative to the 
dimension of the full TO-scheme diagonalization, <i/<i max . 
The ground-state converges very fast while the excited 
states converge somewhat slower toward the exact re- 
sults. To obtain fast converging results for the excited 
states, we directly solve the eigenvalue problem (0} for 
excited energies. The full data points in Fig. [21 show the 
results of such calculations for the first and second ex- 
cited state, respectively. 

Figure shows the squared overlaps of the factoriza- 
tion results with the exact solution versus the dimension 
d of the eigenvalue problem. Excellent results are ob- 
tained for the ground state (which was directly targeted) 
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FIG. 2: Low-energy spectrum for 24 Mg (USD interaction) 
versus the relative dimension of the eigenvalue problem. The 
dashed lines are the exact results, the open data points result 
from targeting the ground-state, and the filled data points 
result from directly targeting the I s * and 2 nd excited state, 
respectively. 

and for the excited states. The inset of Fig.[31shows that 
the angular momenta of the low-lying states are accu- 
rately reproduced even at severe truncations. This indi- 
cates that transition matrix elements can also be calcu- 
lated very accurately. The results of Fig. [21 and Fig. |3| 
clearly demonstrate the efficiency and accuracy of the 
factorization. 

For p/-shell nuclei we use the KB3 interaction 0] and 
compute the ground-states of 56 Ni and 48 Cr. The re- 
spective m-scheme dimensions are c? max = 1.087 x 10 9 
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FIG. 3: Squared wave function overlaps with the exact re- 
sults for the low-lying states in 24 Mg (USD interaction) versus 
the dimension of the eigenvalue problem. Inset: Correspond- 
ing angular momentum expectation values. The results are 
obtained from targeting the ground state. 
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FIG. 4: Ground-state energy E versus the dimension d of 
the eigenvalue problem relative to the m-scheme dimension 
d max = 1.087 x 10 9 for 56 Ni. Inset: Similar plot for 48 Cr 
(d max = 1.96 x 10 6 ). 



and dmax = 1-96 x 10 6 . Figure ^ shows the ground- 
state energies for 56 Ni and 48 Cr (inset) plotted versus 
the relative dimension e?/d max of the respective eigen- 
value problem. We also show an exponential fit of the 
form E(d/d max ) = E + frexp (— cd/<i max ) to the right- 
most seven and six data points, respectively. For 56 Ni 
( 48 Cr) the resulting ground-state energy is only 100 keV 
(30keV) above the exact result while the dimension of 
the eigenvalue problem is dramatically reduced by a fac- 
tor 1000 (factor 12). Figure [5] shows the excitation 
spectrum relative to the ground-state. Our method re- 
produces the level spacings of the two lowest excitations 
very accurately even at the most severe truncation, while 
the spacings to the higher levels are about 300 keV too 
large. Quantum numbers are reproduced accurately for 
d/d max > 0.05. Considering the modest size of the eigen- 
value problem we solved, these are very good results. 
We also compared the m-scheme factorization with a 
particle-hole calculation. The results are listed in Ta- 
ble [n] and clearly demonstrate the fast convergence of 
the factorization method. 

Let us also compare the ground-state factorization pre- 
sented in this work to three related methods. The DMRG 
also bases its truncation on the singular values [l^. Re- 
cent applications to sim ple nuclear structure problems 
were very successful [Til Il5j|. but the application to re- 
alistic nuclear structure problems seems more difficult as 
the convergence is very slow |l6[ . Recently, Andreozzi 
and Porrini |2lj approximated the shell-model ground- 
state by products of eigenstates of the proton-proton 
and neutron-neutron Hamiltonian. This method yields 
quite good results for energy spacings of low-lying shell- 
model states. However, the generation of the eigenstates 
might become impractical for the large proton and neu- 
tron spaces. A third related method is the Exponen- 
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FIG. 5: Level spacings of the lowest excitations with respect 
to the ground-state of the p/-shell nucleus 48 Cr (KB3 interac- 
tion) plotted versus the dimension d of the eigenvalue problem 
relative to the m-scheme dimension d max . The dashed lines 
are the exact results. 



tial Convergence Method (ECM) [ll El [H El . A di- 
rect comparison is not easy since the FPD6 interaction is 
used for p/-shell nuclei, and since ECM results are plot- 
ted versus JT-coupled dimension of the truncated space. 
However, it seems that the convergence of our method is 



considerably faster. For 48 Cr, for instance, our rate of ex- 
ponential convergence is c f=s —31.38 (See Fig.0J), which 
is about a factor eight larger than what is reported for 
the ECM in Fig. 1 of Ref.|12|. For 56 Ni, our exponential 
rate is about a factor 200 larger than the ECM rate [22( , 
and our identification of the exponential region requires a 
m-scheme dimension d w 10 6 (See Fig.^J while the ECM 
requires an m-scheme dimension of 4-5 million |22j . 

In summary, we proposed a new method that factor- 
izes ground-states of realistic nuclear structure Hamilto- 
nians. The optimal factors are derived from a variational 
principle and are the solution of rather low-dimensional 
eigenvalue problems. The approximated states and en- 
ergies converge exponentially quickly as more factors are 
included, and quantum numbers are accurately repro- 
duced. Computations for sd-shell and p/-shell nuclei 
show that highly accurate approximations may result 
from eigenvalue problems whose dimensions are reduced 
by orders of magnitude. 
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